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Abstract 


The pd 7r~^t process has been calculated in the energy region around the 
A-resonance with elementary production/absorption mechanisms involving 
one and two nucleons. The isobar degrees of freedom have been explicitly 
included in the two-nucleon mechanism via vr- and p-exchange diagrams. No 
free parameters have been employed in the analysis since all the parameters 
have been fixed in previous studies on the simpler pp Tr~^d process. The 
treatment of the few-nucleon dynamics entailed a Faddeev-based calculation 
of the reaction, with continuum calculations for the initial p-d state and 
accurate solutions of the three-nucleon bound-state equation. The integral 
cross-section was found to be quite sensitive to the NN interaction employed 
while the angular dependence showed less sensitivity. Approximately a 4% 
effect was found for the one-body mechanism, for the three-nucleon dynamics 
in the p-d channel, and for the inclusion of a large, possibly converged, number 
of three-body partial states, indicating that these different aspects are of 

comparable importance in the calculation of the spin-averaged observables. 
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I. INTRODUCTION 


Pion production/absorption on nuclear systems represents a complex, challenging prob¬ 
lem, and this fact has been known since many years. In the attempt to explain the multitude 
of experimental results collected over more than 40 years by now, many different theoretical 
approaches have been proposed so far, with the aim to improve our understanding of these 
reactions and, more generally, of the hadronic phenomena. 

The simplest approaches which have been employed assume a one-nucleon mechanism 
originated by the vrNN vertex (generally - but not always - a Galilei invariant non relativistic 
reduction of the usual 75 pion-nucleon interaction) and recast the transition amplitude in 
a DWIA formalism where the many-body aspects can be treated with different levels of 
approximations |I|]. First-order corrections from the vr-nucleus multiple scattering series 
leads to the two-nucleon mechanisms where a pion emitted from one nucleon scatters from 
another before leaving the whole nucleus. 

The non-perturbative character of the hadronic interactions, together with the energy- 
momentum mismatch which forbids absorption on single, free nucleons and suppresses one- 
nucleon absorption on bound nucleons, makes these rescattering effects an important aspect 
which cannot be ignored for the reproduction of low-energy data. At the A-resonance, 
however, most of the assumptions used in the early calculations were not applicable, and 
other phenomenological approaches have been developed. 

A relevant success was encountered by the “deuteron” model, originally employed by 
Ruderman [^, where the pd —> yr+t cross section is directly linked to the pp yr+d experi¬ 
mental cross section, by means of suitable nuclear structure functions. This method has the 
considerable advantage to disentangle the question of the choice of the correct interaction 
mechanisms from the knowledge of the bound state wavefunctions. The theoretical uncer¬ 
tainties about the interaction mechanisms at the resonance are simply bypassed through 
the employment of the experimental pp yr+d cross section. The approach has been later 
rehned by several authors, which soon addressed their attention to the main issues such 


2 


as the role of the distortion effects, sensitivity of the results with respect to the available 
model wave-functions, and generalization of the formalism from three-baryon systems to the 
many-nucleon + 1 case. A partial but representative sample of works developed 

along these lines is given in Refs. 

More recently, the analysis of the A{p,'n^)A + 1 reaction in terms of the pp —>■ Ti^d 
process has been extended to spin observables by Falk p, by using the pp —>• amplitude 

analysis of Bugg et al. 0, rather than the averaged cross-section data of earlier works. 
In most cases, an overall fair, qualitative agreement with the trend of the large bulk of 
experimental data has been found, indicating that the deuteron model even in its spin- 
dependent version may be considered a starting point for phenomenological studies which 
include spin-polarization data. This kind of approaches, however, suffers from a number 
of problems which limits possible future rehnements and demands for more appropriate 
theoretical formulations. First, the model considers NN production mechanisms limited to 
the (dominant) 1 —> 0 isospin transition, while a complete model should take into account 
also the effects in the remaining isospin channels. Second, the amplitudes used for input are 
on-shell quantities, while in the theory the off-shell effects should be taken into account. 
Third, there is a certain arbitrariness in the kinematical constraints which dehne the energy 
parameter for the evaluation of the input 2N production amplitude; because of that, the 
results are not uniquely determined. And hnally, in practically all applications the distortion 
contributions are in danger of double countings; indeed, in the evaluation of the distortion 
effects one should subtract the distortion contributions on the active pair, which are already 
included in the pp —>• 7r+d data. 

To overcome these limitations, and to challenge the validity of the “deuteron” approaches, 
there has been a number of papers where more microscopic models have been suggested 


.0|. Here, the dynamical input was not mediated or hidden by the cross-section data, but 


consists of non-relativistic interactions among pions, nucleons, and isobars. In particular, 
the coupled AN dynamics is obtained through the solution in R-space of the two-baryon 
Schrodinger equation. First calculations [§] employed simple s-wave deuteron and triton 
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wave-functions of uncorrelated gaussian forms, and were extremely limited in the number 


of intermediate states included in the calculation. Subsequent analyses P,p!0| used better 
parameterizations (including d states) for the nuclear wave-functions, and increased slightly 
(from 1 to 4 or 5) the number of two-body angular-decomposed intermediate states, hnding 
that these aspects improve the calculated angular distributions without, however, hnding 
an agreement with theory and experiments, especially for the results at large angles, where 
considerable discrepancies persist. One of the conclusions from these studies was that these 
microscopic calculations needed important improvements in order to give a good reproduc¬ 
tion in the normalization and in the angular distributions, and that these improvements 
would lead to enormous complications in the theoretical evaluation of the observables. 

In an attempt to go beyond the description of meson production in terms of single¬ 
nucleon and two-nucleon mechanisms [^, it was shown that meson double rescattering 
could be a good candidate to account for the discrepancies at backward angles, particularly in 
the energy region above the A-resonance. However, the bulk results, while moving towards 
the right direction, were still far from being optimal and this again poses the question of the 
need to overcome the various technical approximations which had to be assumed to keep the 
calculation tractable. The study has the merit to address the attention to the problem of 
three-body mechanisms in meson production/absorption, and this is still an open question 
which we hope will be theoretically disentangled in the near future with the help of the results 
collected in recent years by pionabsorption experiments with large angular detector systems 
|T^^T4[ and by making comparative analyses in the A region with the phenomenology of 
^He photodisintegration |T^ . 


In this paper, we have calculated the excitation function and differential cross-section 
for the pd —> nt reaction in the A region by using single-nucleon and two-nucleon mecha¬ 
nisms. In particular, the two-nucleon mechanism explicitly refers to the intermediate isobar 
excitation with 7r-|-p exchange AN-NN transitions with tensor components. The one-body 
mechanism is mediated by the ttNN vertex, while the two-body process is triggered by the 
ttAN one. The two mechanisms have been decomposed in complete three-nucleon partial 
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waves, while the asymptotic pion-nucleus plane wave has been kept three-dimensional. This 
choice leads to matrix elements with a large number of couplings between different three- 
body states and the technical complications involved has been kept under control with a 
great deal of numerical analysis. Also the intermediate three-baryon (ANN) state has been 
represented in partial waves with inclusion of angular momenta up to / = 2 for the AN 
subsystem. However, in preliminary test calculations H, the quality of the convergence 
has been checked up to / = 5. The antisymmetrization prescription for the three nucleons 
has been fully taken into account via the formalism of the permutation operators and 
brings into the theory further mechanisms which differ from the leading ones by exchange 
diagrams. The matrix elements of the symmetrization permutator P involve a large num¬ 
ber of couplings between different three-nucleon partial waves, and rather than invoking a 
drastical reduction of the number of states, we have accepted to deal with the numerical 
difficulties implied by the approach. 

Following the formalism developed in Ref. |]^ we have embedded these absorption oper¬ 
ators in a Faddeev-based treatment of the three-nucleon dynamics. Accurate bound-state 
wavefunctions have been obtained with high-rank Faddeev-AGS calculations, and a suitable 
generalization of the quasiparticle method to absorption processes has been formulated in 
order to take into account the three-nucleon dynamics in the pd channel. This continuum 
calculation has been performed in the rank-one (separable) approximation and therefore is 
not as accurate as the bound-state wavefunction, but this is the first calculation which, to 
our knowledge, includes a Faddeev-based treatment of the initial-state three-body dynamics 
for the pd ^ 71 + t process (and the other processes related by isospin symmetry). 

Because of the above mentioned reasons, we think that the treatment illustrated in this 
work achieves several improvements with respect to the previous, pioneering, microscopic 
analyses. We have devoted section II to further discussions on these and other important 
aspects of our calculations, such as the treatment of the A mass and width in the inter¬ 
mediate three-body (ANN) Green’s function, and the inclusion of off-shell effects (for this 
second aspect, see also the discussion at the beginning of Sect. III). The remaining part of 
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Sect. Ill compares the results of our analysis with experiments and Sect. IV contains a brief 
summary and the conclusions. 


II. THEORY 


We express the transition amplitude for the ^ pd reaction by the following matrix 
element 


=5{'i'hi'4i'i'Bs)nPo'). 


( 1 ) 


where I'hBs)^ represents the three-nucleon bound-state (BS) and the three-body 

continuum state with ingoing boundary condition and with the deuteron d in the asymptotic 
channel. Both states are assumed to be properly antisymmetrized in the nucleonic coordi¬ 
nates and the state |Pq) is the pion-nucleus three-dimensional plane wave. In a previous 
exploratory calculation [j^ the bound-state wavefunction calculated in Ref. [^] was used. 
That wavefunction originally was determined in Ref. [Q by solving the Faddeev equation 
with the Paris potential and with explicit allowance of the A degrees of freedom and con¬ 
sisted of 48 three-baryon partial waves, of which 22 refer to purely nucleonic states. Because 
of the smallness of the remaining isobar states, only the nucleonic states were retained in the 


calculation of Ref. jig. In the present calculation, we have independently recalculated the 
three-nucleon bound state for various realistic nuclear potentials, and the details of these 
calculations are illustrated in subsection [11 B| . 

The absorption mechanisms are specihed by the operator A. To avoid double countings, 
we avoided purely nucleonic intermediate states in A. This is because we calculate the 
hnal-state interactions among the three nucleons using realistic NN potentials, without 
performing any kind of subtraction in the nuclear potential. 

In this study we consider absorption mechanisms generated by pion-nucleon interactions 
in p waves as described by the ttNA and ttNN vertices. The nonrelativistic interaction 
hamiltonians are 
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HttNA 


( 2 ) 


= /<irp(r)(S-V)(4(r).f), 

niTT J 

for the ttNA interaction, and 

H^nn = f dr p{r){a ■V){${r) ■ t) , (3) 

J 

for the ttNN one. Here, the baryonic density is denoted by p(r), while *h(r) is the pionic 
isovector held. The quantities M and m-,,^ are the nucleon and pion masses, respectively, 
while a and r are the nucleon spin and isospin operators, and S and T are the corresponding 
generalization to the isobar-nucleon transition. 

The ttNN vertex dehnes the simplest one-body absorption mechanism and is sometimes 
referred to as the impulse approximation (lA) mechanism. However, this is suppressed be¬ 
cause of energy-momentum mismatch, therefore two-nucleon mechanisms dominate. These 
are taken into account through the A-rescattering process, where the ttNN-NN inelastici¬ 
ties are modeled through a ttNA-AN-NN two-step transition. Here the hrst transition is 
triggered by the ttNA vertex given above, followed by an intermediate AN propagation and 
by the AN-NN transition. The intermediate propagation of the isobar is described by the 
Green’s function 


Gn — 


E + M - Ma - 


( 4 ) 


2/ia 2i^A 

where E is given by the pion energy plus the target kinetic energy in the c.m. system. 
The Jacobi variables p and q are the pair AN and spectator nucleon momenta, respectively. 
Correspondingly, pA and ua are the reduced masses of the pair and spectator-pair systems. 
AJa ~ Af is the isobar-nucleon mass gap. Since the A-isobar is not a stable particle but a 
resonance, its mass is endowed with an imaginary part Ta, associated to the decay width of 
the isobar 


MA = MA-Es-'-TAiE), (5) 

with Ma = 1232 (MeV) and Eg being the energy-shift parameter. The energy dependence 
for T has been modeled phenomenologically in a previous analysis of the Trd <->■ pp reaction 
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T^. Here we use the same analytical expression which has been obtained in that analysis, 
with the Fa energy dependence extracted by the condition 


a(E)^ 


D 


( 6 ) 


(E - Bk)2 + ’ 

where <j{E) is the experimental excitation function for pion absorption on deuterons. In the 
numerical calculation, the £t to cr{E) obtained by Ritchie has been used. 

In a meson-exchange framework the NA-NN transition potential can be ob¬ 

tained from the tt- and p-exchange diagrams 


VwA = (VJi + • T,) 

( 51 * • • Q)GM 

The operator Q is the baryon-baryon transferred momentum, nip is the mass of the p- 
meson and the function ^(a;), which describes the meson-exchange propagation, takes into 
account the mass difference between nucleons and isobars 


= TT-y + 


2uj‘^ — M + (jj) 

1 1 


( 8 ) 


2a;2 2a;2 + - M) ' 

Here, uj is the relativistic energy of the exchanged meson and, as shown in the last expression, 
we have taken into account the AN mass difference in an approximated way in order to 
obtain analytical expressions when performing the partial-wave expansion. Each meson- 
baryon coupling in the transition potential has been endowed with a phenomenological form 
factor of monopole type, with the exclusion of the pNA coupling, where a dipole-type form 
factor was assumed. 

Finally, we discuss how the three-nucleon dynamics can be incorporated into the theory, 
or equivalently, how we calculate ^|. The method is similar to previous procedures 

developed to incorporate hnal-state interactions (ESI) in the photodisintegration of three- 


body systems P7| ,^ |23H^. First, we introduce the operator P 23 which exchanges spin. 
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isospin and position coordinates of nucleons 2 and 3. We then introduce the cyclic and anti- 
cyclic permutation operators P 2 and P 3 respectively. They exchange the global coordinates 
of the three nucleons in the following manner 123 —231 and 123 —312, and can be 
expressed in term of the pair exchange operators as P 2 = P 12 P 23 and P 3 = ^ 13 ^ 23 - This 
leads to the full permutator 


P = P2 + P3, (9) 

and the normalized symmetrizer S = It is now possible to derive the complete wave- 
function in terms of the asymptotic channel wavefunction (the subscript “1” denotes 

antisymmetrization with respect to the pair made of nucleons labeled “2” and “3”) 

s{¥,-^\ = ,{^d\{l + TGo)S, ( 10 ) 

where the three-body operator T satisfies a Faddeev-like equation 

T = Pt + PtGoT (11) 

{t is the two-nucleon t-matrix). 

We now dehne the Mpller operator = 1 -|- TGq, in which case 

= ( 12 ) 

This operator satishes the Faddeev-like equation 

= 1 + , (13) 

and by its means the pion-disintegration amplitude can be rewritten as 

2l, = i($,|fl(-)t,S^|vI/Bs,Po)- (14) 

Because of the two equations above, the full amplitude satishes the integral equation 

A, = i($,|F^|Tbs, Po) + , Pq) • (15) 
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The first term on the right-hand side corresponds to the plane-wave contributions, while in 
the second term the NN rescatterings (to all orders) are singled out. 


As is well known, when the two-body transition matrix is represented in a separable form, 
the Alt-Grassberger-Sandhas (AGS) equations for neutron-deuteron scattering reduce 
to effective two-body Lippmann-Schwinger equations. The same happens for the calculation 
of hnal-state interaction effects in the photodisintegration of the triton (7 + f —>■ n + d) 
p7|,|^,^3|-^ , where one obtains a similar effective two-body Lippmann-Schwinger equation, 
the only difference being that the driving term of the n-d scattering equation (he. the 
particle-exchange diagram, the so-called “Z”-diagram) is replaced by the off-shell extension 
of the plane-wave photoabsorption amplitude. Here, we use a similar scheme for pion 
absorption on the three-nucleon systems. 

In order to accomplish this, we use the separable expansion method proposed by Ernst, 
Shakin, and Thaler (EST) for representing a given NN interaction. The EST method is 
very reliable and has been tested in the past for bound-state and scattering calcula¬ 

tions P5|-|H7]. In this approximation the original interaction is expressed in separable terms 
of the form 


N 




(16) 


where N is the rank of the approximation. The detailed structure of the potential, the 
parameters for the formfactors |/^) and the interaction strength A can be found in Ref. 
0 , |37[| . The input two-body transition matrix t is given as a series of separable terms 


t = E l/.> A.- (W ■ (17) 

flU 

To simplify the notation, we restrict the sum above to just one value for the indices fi 
and u, so that the transition matrix takes the rank-one form t = |/i) A (/i|. In A as well 
as in the form factor |/i), we have also omitted the proper energy dependence, but it must 
be remembered that when the two-body f-matrix is embedded in the three-particle space 
the correct energy dependence is upon E — (||j), where the energy of the spectator nucleon 
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has to be subtracted. The separable representation reproduces the correct negative-energy 
bound-state pole of the two-body t-matrix (for the deuteron quantum numbers) if the form 
factor satishes the homogeneous equation 


VG„(Ei)\S,) = |/i> 


(18) 


and, within a normalization factor, Go{E(i)\fi) is the deuteron wave |$d). 

Use of the separable representation of the t-matrix input in the integral equation (|T^ 
leads to the effective two-body equation 


A{q, E) = B(q,E)+I q'^dq' V{q, q', E) A{E - A(q', E) 


(19) 


with the dehnitions 


A{q,E) = i{/,,,|Gi+>li>->tS.4|$Bs,P;), 
B(<(.S) = i(/.,<(|GhS.4|>l<Bs.P;). 

V{q,q',E)^l{f„q\G^;>P\fl,q')l. 


( 20 ) 


Here, A and B represent the off-shell extension for the full and plane-wave pionabsorption 
amplitudes, respectively, while V is the effective two-body potential which represents the 
one-particle exchange diagram between different sub-cluster rearrangements. 

These amplitudes are decomposed in three-nucleon partial waves, while the pion-nucleus 
incident wave is treated in three dimensions. Details on the employed representation were 
given elsewhere [^]. The representation of the three-body states is dehned in momentum 
space and the partial-wave decomposition is discussed within the jl coupling scheme. The 
index a refers to the whole set of quantum numbers (he. orbital momentum, spin, total 
angular momentum, and isospin of the pair {ls)j; t, of the spectator (Act)/; r, total angular 
momentum, isospin and associated third components JJ^]TT^). 

The absorption mechanisms we include in our calculation have the following structure: 


B^{q\a\E) = 2 V 3 
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and 


X ^ / p'‘^dp' q"‘^dq" p\dp/^p^dpq^dq 


X 


cx." 

{p''q''a''\VN/^\PAq''a/\) 


h{a\p') {p'q'a'\P\p"q"a") 

E - 32 ^! _ 2 ^ + 

^ AM M ^ 


{pAq”aA\H^NA\pqa) {pqa\^Bs) , 


B^{q',a\E) = 2V3 
X H / P^dp pidpAP^dpq^dq - - 

OL/^,OL'^ ^ -— — 

{p'q'a'\VNA\pAq'oLA) 


E -^al + ie 

^ AM M ' 


{M-M, 


q'2 Pa 


^ {PAqaA\H^NA\pqoi) (pgald'Bs) 

'a 


2i/a 2/iA 

for the A-rescattering mechanisms, and 

Bi,^iq\a',E) = V3 
X J p‘^ dp p"'^ dp” q'^dqp^dpq^dq 

Ol" ,0L 

X {pqa\H^MN\pqoi) {pqa\^Bs) , 


h{a\p') {p'q'a'\P\p”q”a”) 

E - lal - sl + ie 
^ AM M ^ 


( 21 ) 


( 22 ) 


(23) 


B’a*i^,a',E)^V3 

X ^ / p^dp p^ dp q^dq 


fi{a',p') {p'q'a'\H^NN\pqa) {pqa\^Bs) 

E - ^a!l - e! 1 A_j(: 

^ AM M ' 


(24) 


for the one-body lA mechanisms. Clearly, all these mechanisms add coherently to give the 
total amplitnde. The snbscripts “E” and “D” differentiate between exchange and direct 
mechanisms, respectively. The “D” mechanisms correspond to the probability that nncleon 
“1” is the free nncleon in the ontgoing channel (hence nncleons “2” and “3” form a denteron) 
while the other two cases (he. nncleon “2” or “3” as free ontgoing particle) are assembled 
together in the “E” mechanisms. These exchanges in the rearrangement channel are per¬ 
formed by the permntation operator P. We refer to the appendix of Ref. for details 
on the partial-wave matrix elements for this operator. Each A-rescattering mechanism is 
mnltiplied by the factor 2a/3 which arises from the mnltiplicity of the possible two-body 
mechanisms and by the normalization factor dne to Panli principle. The same considerations 
lead to the factor a/3 in the impnlse approximation. We also observe that a conpled-channel 
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structure has to be intended for the deuteron form factor since the deuteron d- 

wave component is obviously taken into account in our calculation. However, these coupled 
orbital-momentum components of the pair must be summed up coherently. 


We assume that such a sum (over the coupled Vs) is performed at the present stage, so 
that the index a', from now on, collectively denotes the set of quantum numbers s, j, t, A, 
I, r, J, J^, T, and T^. 

Finally, we observe that the direct contribution to the A-rescattering term, Eq.(^), 
vanishes on-shell. (On-shell, the two-nucleon quantum numbers s, j, and t are hxed by 
the deuteron state, and the relative motion of the spectator q is hxed by total energy 
conservation.) This vanishing effect is due to isospin considerations, since the intermediate 
AN pair cannot be directly coupled to a hnal deuteron-like (t = 0) pair. It does couple, 
however, through the permutation operator P in the exchange contribution, as well as 
through hnal-state interactions, where one can well have isovector pairs in the intermediate 
states. 

Once the total absorption amplitude E) has been obtained, the pion- 

absorption excitation function is given by 



(25) 


with phase-space factor 



(26) 


where 





Eii+pf 


E„ = ^/WT^ 


( 27 ) 



E^°^ — E]\f + Efi — Et^ + Et . 
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Here, Md, Mt, are the deuteron and three-nucleon masses, respectively. The momentum q 
is the on-shell momentum (in c.m. frame) of the outgoing nucleon. 

Similarly, the unpolarized differential cross-section for the 7r~^t pd reaction is given by 


^ y X 
dil ^ 2 


C(A/iV;A‘A'V‘)rf(fi) 

J,J^,X I 


1 - K 

2 J, (28) 

X J I 


where K is the channel-spin of the deuteron-nucleon system, i.e. the sum of the spins of 
the two fragments K = j + a. 

Finally, detailed balance gives the expression for the unpolarized cross section for the 
inverse (pion production) reaction, 

“I absorption 

. (29) 


dc , 

production - / \ 2 

da , , 





where 9 is the c.m. angle of the deflected particle. 


A. Partial—wave absorption amplitudes 

We outline here the basic ingredients we used for the calculation of the absorption ampli¬ 
tudes in partial waves. The section is mainly technical and can be ignored on a first reading, 
if one is not interested in the details of the calculation. 

For brevity, out of the four mechanisms discussed previously, we have selected only 
two mechanisms, namely the exchange A-rescattering, which is the dominant one, and the 
direct LA. term. The remaining mechanisms, direct A-rescattering and exchange LA. have 
a similar structure and do not introduce any novelty. 

To avoid unnecessary complications in the formulas, we give here the amplitudes on shell, 
i.e., we have substituted |<hrf) in place of Go{Ed)\fi) and have denoted by uisj{p) = u{p) the 
s and d component of the deuteron in momentum space. For the inclusion of the three- 
nucleon dynamics via the AGS equations, these same amplitudes have been extended off 
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shell (see Eq. (p!8|) ). To fix the notation, the exchange A-rescattering amplitude in partial 


waves can be denoted as follows 


^ 2a/3 ,q,PS) i 


(30) 


where q' is the c.m. momentum of the N-d system, Pq is the c.m. momentum of the 
t channel and denotes the beam axis or equivalently the z-axis. While the former is a 
one-dimensional variable in a fully decomposed partial-wave scheme, the latter is three- 
dimensional, since we decompose in partial waves only the baryonic coordinates, not the 
pionic one. The operator Ai represents the ttA-A transition on nucleon labeled “2”, in¬ 
termediate A propagation, and a AN-NN transition between nucleons “2” and “3”, while 
nucleon “1” acts as spectator. The process is diagrammatically illustrated in Figure 1 of 


Ref. IT 3 , while P is the three-body permutator previously introduced. For the detailed 
form of P in partial waves, we used the expression given in the appendix of Ref. . 


With a' we collectively denote the quantum numbers for the jl decomposition of the p-d 
channel listed according to Eq. (3.2) of ref. [^, while with a the quantum numbers for the 
three-nucleon bound state in the LS scheme are assumed. The details of the calculation in 
this scheme can be found in Ref. []^|, while similar calculations in the jl scheme are found 
e.g. in Refs. |^,|M[|. 

The resulting expression is 


{u, a, q'lPAil'^BS, a, Pq) = 


E 

L/\ S/\ Sz'P^ 


g{a'a"khy[l'2) T{a"a) 


X B{a"a/ 2 ^aLi^Si^Sz'm) Xi{oi a'a^aklil'iL^S^Szm ), 


(31) 


where 


^ ^ iPpiSM + 

727T^ml{M + 


( 32 ) 


~g {aa'khhl'iQ = -isjfxil'sj't'X'I' 


1 

2 

T 


t 

t' 
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oo 1 oo 

li{a'a'aAakkl'iLASASzm) = jp'^dp'uvs>j'{p') J dx J, Jp^dp 

-1 ^ ^ 0 


0 


1 


X / dcosP 0 ;^^(cosPa) 0rm(cosP) 

^ _ Q _ 

2/ia 2i/a 


-1 


X / dcos Q' 0A"n(cOsQ') 0 a„(cOsQ) d'< 5 (p, g) . 


-1 


In this last equation, all the relevant momenta are defined in terms of the channel 
menta q' and Pq and of the integration variables through the formulas 


(33) 


(34) 


(35) 


(36) 

mo- 
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p" = \/\p’‘^ + + Ip'q'x^ 


€ = ^p'"^ + - p^q'x : 


/, 3M + ^ „ 3M + ^ ' 

Pa = \ ijrvj —;;- Po)^ + P^ + 2-77—z- Pq pcosP 


cos Pa = 


'6M + 3m-n 
m^P^ + pcosP 
PA 


6M + 3m^ 


cos Q = 


q = \J + 9*^ + \Po Q*^ cos Q', 
^Pq + q* cos Q' 


(37) 


In Eq. (^) Pk{x) represents the Legendre polynomials, 0;m(cos6*) denotes the asso¬ 
ciated Legendre fnnctions (normalized according to Ref. 0), Qli' {p,p') represents linear 
combinations of Il-kind Legendre polynomials (see Eq. (43) in Ref. |]^), originated by the 
7r-|-p meson-exchange diagrams in the tensor AN-NN force, and hnally 4/a(p, q) is the triton 
wavefunction in momentum space and in the LS scheme. In the same equation, it is also 
assumed that n = — m, while the identities h + h = I' and l[ + 12 = X' hold in both 

Eqs. (^^ and (|36|) . 


Similarly, the direct LA. amplitude in partial waves can be expressed as follows 


Bi"'- = i/3{uaV|B,A,K|1'Bs.a,PS>. 


(38) 


with the ttNN vertex acting on nucleon “2” as a one-body operator. 
The calculation of this matrix element yields 


{u, a , q\H.„NN\^BS, «, Pq) = P2 


B 2 {aaL'S'S,m) X 2 {a ah'S'S,m ), 


nJr^'^TlOl 

L'S'S^m 


(39) 


where 


P2 = 


iP^{3M + m^)f^NN 


(40) 


B 2 {a'aL'S'S,m) = (-)i+®+^ +*+^ +5 I s' f t' S T j' I' L' S' 

] 


i i t' 


1 i i 
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X c[l'X'L'- m{J', -S,- m)(j; - S,)) C{r\L- m{J'^ -S,- m){J'^ - S,)) 


(41) 


and 

OO 1 

X2{a'aL'S'Szm) = Jp'^dp uiig’j'ip) J dcosP'0z'm(cosP') ©(-^(cosP) 

0 -1 

1 

X J dcosQ'Qx'n{cosQ')Q-xn{cosQ)'^a{p,q) ■ (42) 

-1 

As in the previous case, all the relevant momenta have to be expressed in terms of the 
channel momenta q' and Pq and of the integration variables 



, 3M + „ 3M + m,r , vs,, 

P + -Pq") - 2 ^,. , cosP' 


6M + Sm-!, 


QM + StUt, 


q'^ + (1^0 f + IPq q'^ cos Q', 
^g^Po-+p'cosP' 

p 


cos Q = 


\Pq + q' cos Q' 


q 


(43) 


B. 3N Bound—state calculation 


For the calculation of the three-body bound-state energy Et and wavefunction |\hBs) we 
start from the Lippmann-Schwinger equation 

|4/bs) = G^(-^t + *c) |4 /bs) • (44) 

Using the resolvent equation 


G{z)=G0{z) + Gp{z)V0G{z) 


and performing the e-limit we end up with 


(45) 
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'I'bs) — Gi3{ET)(y — Vg)|^'Bs) 


(46) 


where V is the total interaction snmmed over the pairs, Vp is the nnclear interaction of the 
pair labeled by /?, Vfs denotes V — Vp, and Gp is the channel resolvent for the two-cluster 
partition. 

If we now introduce the position \Fp) = {V — Vp) Ift'es) and the relation VyGj = TjGq, 
we obtain the celebrated equation 

|Fs> = i:(l-is,)r,G„|f,). (47) 

7 

Here, the summation runs over all two-fragment partitions 7 . The “form-factors” \Fp) are 
related to |4 /bs) by 

I^Bs) = E G'o T, Go \F,) = Y. 1 ^ 7 ) > (48) 

7 7 

where the are the standard Faddeev components. This can be shown by using the 
dehnition of \Fp) and again the relation V^G^ = T^Gq'- 

\^^) = GoT^GolF^) = GoV^G^lF^) = GoV-,G^{V -V^)\^!bs) = GoV^l^jip). (49) 

The bound-state equation (|47|) has to be projected onto the three-body partial waves 
\pqlhJT). In this subsection we use the so called channel-spin coupling. The label h denotes 
the set of quantum numbers {'qKX) , where K and A are the channel spin of the three nucle¬ 
ons [ with the coupling sequence {j, a)K ] and the relative angular momentum between the 
two-body subsystem and the third particle, respectively. In this paper, the channel-spin 
coupling K has been already introduced in Eq. (PB]). Here, the total angular momentum of 
the three-nucleon system J is given by the coupling sequence (X, A) while T has been al¬ 
ready dehned as the total isospin. I is the usual orbital angular momentum, and rj = [s, j]t) 
collectively denotes the spin, angular momentum [ with the coupling sequence {l,s)j], and 
isospin of the two-body subsystem. Dehning 

00 

= E / dPP" Mp) {pqlbJT\Go\Fp) , (50) 

I n 
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inserting the finite-rank form for the T-matrix, and using the Pauli principle, we obtain the 
one-dimensional equation 


OO 

= E Jdq'rf-AV'^'iq.q'. Et) (Et - F”' (q') , 


( 61 ) 


with 


OO OO 

^V“'(g, g', i?T) = 2 ^ / dpp^ dp'p'^ flip) {pqlbJT\Go\p^fn^JT) /,(p') . (52) 

0 0 

The recoupling coefficients for the coupling scheme herein employed can be directly found 
in Ref. [Q , or alternatively, they can be obtained starting from other coupling schemes 
by means of the usual transformation algebra. 

Equation (|5TD can be treated as an eigenvalue problem. The energy is varied until the 
corresponding eigenvalue is 1. The ranks for the partial waves taken into account for the 
bound state calculations for the different potentials are listed in Table B Results for the 
binding energies can be found in Table B 

The whole wavefunction can now be calculated by either using Equation (^8|) or by 
applying the permutation operator P defined in Eq. (P) on one Faddeev component |2^ 


|Tbs) = (1 + P)|^i). 


(53) 


This second representation has the advantage to be explicitely independent from the T- 
matrix representation, and is therefore computationally more convenient when the rank of 
the representation becomes large. The high accuracy of the wavefunction has been shown 


in Refs. 33,34 


For the calculation of the Bornterms as illustrated in the previous section, the wavefunc¬ 
tion was transformed into the LS scheme via 


{iilX)L,isa)S)jr\ 

_ , , , , _ S k] S K 

j2i-iy+^+^+^+^+s+^^jLSK < > < 

jK (j j s J A L 


mis)ja)KX)J,P 


(54) 
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C. 3N continuum calculation 


In this section we briefly summarize the continuum calculation. As already mentioned 
the method used here is similar to the techniques developed to incorporate final-state inter¬ 
actions (FSI) in the photodisintegration of triton . The only difference is the 

inhomogeneity, i.e., the Bornterm, of Eq. (0). Thus, with this replacement our working 
program for photodisintegration has been easily modified for pionabsorption. Both pro¬ 
grams are based on a program for n-d scattering |^, which has been rewritten for the 
EST-method used here and the above replacement. The Faddeev results for the continuum 
presented here are all calculated using a rank-one representation of the two-body t-matrices 
in the kernel of the integral equation. In the case of photodisintegration it has been shown 
3^ that this is enough for the observables calculated. 

In the case of the pionabsorption the Bornterm is much more complicated and has more 
structure than in the case of photodisintegration. To achieve the required accuracy we used 
70 grid points in the q variable for each channel for the off-shell extension of the Bornterm. 
For details of the numerical solution of the integral equation we refer to Ref. . 


III. RESULTS 


Our analysis of the pd vr+f reaction is essentially a parameter free analysis, because all 
the parameters of the model were fixed in previous studies on the simpler vr+d —>■ pp 

reaction. This was possible because the character of the calculation is sufficiently microscopic 
that the tuning parameters are basic quantities such as coupling constants and cut-offs at the 
meson-baryon vertices. Obviously these very same parameters enter in both pp yr+d and 


pd <->■ reactions. We refer to the analyses of Refs. ||T^J^ for a more detailed discussion 
on the TT+d —pp reaction; here we limit ourselves to summarize few aspects which are 
specific of this approach. The coupling constants for the A-rescattering mechanism are the 
ones referring to the nNN, ttAA, pNN and pNA vertices, which were taken from Table 
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B.l of Ref. (model III). The cut-off parameters (in MeV) of the four vertices are 1600, 
900, 1200, and 1350, respectively, which are practically equivalent to the ones given in the 
same Table B.l cited above. 


In the analysis of pion absorption on deuterons which we used as input for our study, there 
were only two authentically phenomenological ingredients. First, the way the A-resonance 
propagates in the intermediate states, and, second, the off-shell nature of the ttNN and 
ttAA vertices when the pion is the asymptotic particle, and hence rigorously standing on 
its mass shell. 

Concerning the hrst point, the resonance peak in the experimental pion-deuteron ab¬ 
sorption cross-section is below the position of a free A (as observed by tvN scattering) by 
a few tenths of MeV. This effect could not be explained within the dynamical model, and 
therefore the A resonance peak has been positioned downwards by introducing an energy 
shift parameter. Phenomenologically, the magnitude of this shift is of the order of 30-35 
MeV; it has, however, to be increased of about 20 MeV, (see Refs. [^,^]) if the motion of 
the other nucleon is taken into account. As for the isobar width, this was fixed to 115 MeV 
at the resonance, and therefore corresponds to the value of the free resonance width. Note, 
however, that the resonance peak for pion absorption on deuteron is broader, being around 
To = 150 MeV This broadening of the width is fully described by the nd —>■ pp 

model calculation, once all the mechanisms, including FSI, are taken into account. 

The second point with a certain deal of phenomenology concerns the off-shell nature of 
the pion-baryon vertices, in case the pion is the asymptotic particle. When the pion is on 
its mass shell, the extended structure of the nNN and vrVA vertices may be governed by 
the nucleon momentum, if the nucleon is the off-shell particle. The importance of including 
these baryonic off-shell effects has been pointed out also in a recent study of nN scattering 
Thus, the two vertices have been endowed with a form factor depending upon the 


momentum k of the interacting nucleon. 




(66) 


N,A 
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For the two-nucleon system corresponds to the square of the Jacobi coordinate p, 
while for the three-nucleon system, an averaging over angles yields ~ + 0.25g^. 

At the level of the two-nucleon system, the model has been tested against the enormous 
variety of experimental data collected for the •«-*• pp reaction, and its strength and 
weakness, especially at the level of spin observables, have been thoroughly discussed in Refs. 

WM- 


To exhibit the quality of the results at the level of the two-nucleon system, we report 
in Figure 1 the total production cross section for the pp reaction from threshold 

to above the A resonance. The parameter p corresponds to the pion momentum (c.m.), in 
units of pion mass. The experimental data were extracted from a collection of Refs. |46|-|^. 
All data have been converted into production data using the principle of detailed balance, 
if necessary, and in case of the data of Ref. |^, charge-independence considerations were 
applied. The full line is the result obtained with the Bonn B potential, and with the p-wave 
pion-nucleon mechanisms which, as discussed in the previous section, include both isobar 
and non-isobar degrees of freedom. Practically indistinguishable from the full line there 
is a dotted line corresponding to a similar calculation but with the Paris potential. The 
differences between the two curves slightly increase with energy and at the resonance peak 
the Paris result is less than 5% smaller than the Bonn B one. 

The other two curves (dashed and dashed dotted, respectively) correspond to the ad¬ 
ditional inclusion of s-wave pion-nucleon interactions in the hnal state of the production 
process (or conversely in the initial state, for pion absorption). The leading contribution for 
this 'kN-'kN interaction was obtained via a p-meson exchange model in Ref. [^. The two 
curves differ in the choice of the NN potential, the Bonn B being referred to by the dashed 
line, the Paris by the dotted-dashed curve. 

The log scale emphasizes the threshold results, and shows that in order to achieve a 
correct reproduction of data over three orders of magnitude, the pion-nucleon interactions 
via p-meson exchange have to be considered as well. Note, however, that the process around 
the A resonance is dominated by the p-wave mechanisms which are triggered by the nNN 
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and tiNA vertices. On a linear scale, the threshold effects in pion prodnction dne to the 
s-wave mechanisms wonld appear smaller. 

At this stage of onr stndy of the pd •«-*• 7r~^t process, we have exclnded the p-mediated 
pion-nncleon interaction in s-wave, and limited the analysis aronnd the A resonance. Fnr- 
ther work is needed to inclnde the low-energy mechanisms, and to extend the calcnlation 
at the prodnction threshold. In what follows we considered pion-nncleon interactions in 
p-waves only. However, both A and non-A degrees of freedom were considered. We note, 
however, that we treat the pion-nnclens wave three-dimensionally and carefnlly calcnlate 
the relevant kinematical transformations from the pion-nnclens to the pion-nncleon vari¬ 
ables. This reqnires additional integrations over angnlar variables in order to prodnce the 
absorption amplitndes (see Eqs. (|36|) and (^^ ). And, becanse of this, the total angular 
momentum J is not restricted, neither is the parity P, and thus we have taken into account 
a few states: J = (|)^, J = (|)^, J = (|)^; and J = (|)^. For each of these three-body 
quantum numbers, one should consider in principle an infinite set of partial waves for the 
two-body subsystem, which matches a corresponding infinite set of partial waves for the 
spectator particle to give a fixed . Clearly, a truncation over a limited number of states 
is necessary. Table 1 shows the smallest set of three-body partial waves which have been 
included in our calculation. Here, four two-nucleon states have been included, namely 

^Pi, and ^D 2 , while the maximum values considered for the orbital and angular mo¬ 
mentum of the spectator were A = 5 and I = ^, respectively, for a total amount of 82 states. 
Such three-body states are indicated in the table according to the notation A/. 

All our calculations, unless otherwise explicitly indicated, have been performed with 
a much larger set. In particular, in calculating the matrix elements of the absorp¬ 
tion/production matrix elements with exchange operator, we took into account 18 two-body 
partial waves ^Pq, ^Pi, ^P,, ^Pa, ^Pg, 2 ^ 3 , ^Pg, 2 ^ 4 , 1 ^ 4 , 

^P 4 , and ^^ 4 . This, with the same cuts in the spectator quantum numbers A and / specified 
as before, yields 34 channels for = (|^), 58 for = (|^), 70 for = (|^), and 70 for 
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= (|^) for a grand total of 464 three-body states. In solving the Faddeev eqnations for 
the three-nncleon dynamics in the p-d channel, of the 18 nncleon-nucleon waves we have 
inclnded the hrst 10 states {j < 2) since for these the NN state-dependent interaction was 
available. 

Fignre 2 shows the integral cross section (in pb) of the pd —> yr+f reaction aronnd 
the resonance. The x-axis dependence is npon the parameter p, which has been previonsly 
dehned. The experimental bars were obtained with the help of a collection of data contained 

for 7r“ 

absorption on ^He at 64 and 119 MeV, assnming time reversal and charge symmetry. The 
full, dotted, and dashed lines differ among each other for the NN potential which has been 
nsed to generate the asymptotic bonnd states, i. e. the denteron and triton. The three cnrves 
represent calcnlations with Bonn B (solid line), Bonn A (dotted), and Paris interactions. 
We hnd that at the peak the Paris resnlts are smaller than the Bonn B resnlts by almost 
25%. This relative difference is snbstantially larger than the 5% difference obtained for the 
simpler pp —> n^d reaction. Experimentally, the pd —> cross-section is smaller by a 

factor 80—85 with respect to the pp —>■ 'w'^d one. This large snppression, dne to the small 
overlap of the denteron wavefnnction in the incoming channel with the pion prodnction 
matrix elements, is fnlly reprodnced by onr calcnlations. 

The nnpolarized differential cross section is reported in Fignre 3 for a variety of energies 
spanning the A resonance. While Fignre 2 shows the normalization of the cross section, in 
Fignre 3 we have addressed onr attention to the pnre angnlar dependence and therefore all 
the cnrves are normalized to the experimental data. Thns, the three cnrves referring to 300 
MeV (lab energy of the proton) have been mnltiplied by a factor. The fnll line, referring to 
the resnlt with Bonn B potential, has been mnltiplied by 1.075, the dotted line (Bonn A) 
by 1.016, and the dashed line (Paris) by 1.329. The three factors are dne to the differences 
in normalization immediately perceived in Fignre 2 for r] = 1, which is the corresponding 
valne for that energy. 


in Refs. 


To this collection, we have added the experimental resnlts of Ref. 


25 




At the resonance peak, the angular differences between the various NN potentials are 
practically zero. These differences, however, increase in moving away from the resonance 
in both directions. We find larger differences above the resonance, for backward angles. In 
particular, calculations with Paris potential seem to reproduce better the data at backward 
angles. 

In the first sector (E = 300 MeV) of the figure, the data were taken from the inverse 
absorption reaction, vr” ^He —> nd at 64 MeV, Ref. [^, using the detailed balance and 
charge symmetry. In the same picture we have reported also four points obtained from Ref. 
0 for the production reaction at 305 MeV. The four points are clearly visible because they 
stand above the rest of the data set. From the same reference 0 , we have taken also the 
data at 330 MeV, shown in the second sector, while the data at 382 MeV were obtained from 
pion absorption reactions at 119 MeV. In particular, data from both reactions 7r~ ^He — nd 
(Ref. [^]) and vr+t — pd (Ref. [^) have been included. The data at 450 and 500 MeV were 
taken from Ref. using charge independence, and finally, at 605 MeV, we considered the 
older data of Ref. [^. In terms of integrated cross-section, that last set of angular data 
corresponds in Figure 2 to the datum aX rj = 2.4, which has the largest error bar. Because 
of this, the large multiplicative factors we found at this energy (5.956 for Bonn B, 5.413 for 
Bonn A, and 7.982 for Paris) seem to be attributable more to normalization problems in 
the data, than to the model calculation. For all the remaining energies, the normalization 
factors were all well around one. 

Table 2 compares for the three NN potentials the integrated cross-section calculated 
under various conditions. The first row exhibits the contribution arising from the sole ttNN 
vertex (denoted lA, Impulse Approximation). In the second row (PWa), we show a plane- 
wave calculation which includes the isobar degrees of freedom via the ttNA vertex. In the 
third row, PWb, the number of intermediate NN states have been increased from 4 up to 
18 two-body partial waves, and from the sole S state up to the D states in the AN orbital 
momentum. The effect in the total cross section is about 3%; however, as shown in the next 
figure, there is a not large but sizeable change in the angular dependence. In a forthcoming 
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article we will show that certain spin observables like AyQ are extremely sensitive to the 
number of waves included in the intermediate states. 

Finally, in the last two rows (ISIa and ISIb) the effect of 3-nucleon dynamics in the 
incoming channel has been taken into account with the Faddeev formalism. Similarly than 
in the previous two rows, these differ between each other for the number of intermediate 
states. In both cases, the effect of the three-nucleon dynamics yields an increase of about 


4% or less. Around the same value for p, the pd —>■ cross-section extracted from Ref. jHO 


gives 35 ± 3.2 pb at p = 1.32, while data extraction from gives 32.7 ± 8.0 at = 1.43. 

Figure 4 shows the differential cross-section calculated for p = 1.36 in various conditions, 
in comparison to the data obtained for tt" production at 350 MeV assuming charge inde¬ 
pendence 1^. The solid line is the solution of the Faddeev equation for the three-nucleon 
dynamics in the incoming channel. The corresponding plane-wave results are shown by the 
dotted line. The comparison shows that the three-nucleon dynamics have a relatively small 
effect on the differential cross-section. On a linear scale, the effect is more pronounced at 
forward angles, however, on a log scale it becomes evident that the overall effect is simply a 
rescaling of the curve, without changing the angular dependence. The dashed line shows the 
results obtained with a limited number of channels (a total of 82 instead of 464). Differences 
are seen at both forward and backward angles. It is therefore important to consider conver¬ 
gence with respect to the number of three-body states included in the calculation. Finally, 
the dashed-dotted line contains the effect of the sole nucleonic intermediate states. For rea¬ 
sons of visibility, these non-isobaric effects have been multiplied by a factor of 10. Overall 
the effect is small, but at backward angles its contribution is larger than 10%. The calcu¬ 
lations shown in the hgure were performed with the Paris interaction. Practically identical 
angular dependences have been obtained with the Bonn B and Bonn A potentials. 

Finally, Figure 5 exhibits the dependence of the differential cross section upon the pa¬ 
rameter p for forward and backward angles. Due to the smallness of the three-body effects, 
the calculations have been performed in plane-wave approximation. The theoretical results 
have been divided by 2 for comparison with the experimental data for vr® production ob- 
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tained at Saclay ||^. The Bonn B and Paris curves have been normalized at E = 350 MeV 


to the data of Ref. 


1601 


IV. SUMMARY AND CONCLUSIONS 

We have analyzed the pd —>■ reaction around the A-resonance with a model calcu¬ 

lation which explicitly includes isobar degrees of freedom and meson-exchange diagrams. 
The elementary production/absorption mechanisms were tested on the simpler pp —>■ n~^d 
reaction. In particular, the position and width of the isobar resonance were modeled to re¬ 
produce the excitation function of the pp —>■ 7r~^d process, and off-shell effects in the baryonic 
coordinates were taken into account in both vrNN and ttAN vertices. The same absorption 
mechanisms (without further changes) have been embedded in a Faddeev-based treatment 
for the pd —>■ process, where the three-nucleon bound state and the three-nucleon 
continuum dynamics in the initial channel have been calculated using the Faddeev-AGS 
formalism. The computational method is similar to the one employed recently for the triton 
photodisintegration. We have checked our results against integral and differential cross- 
section data in the resonance region. Ending that the unpolarized experimental data are 
reproduced reasonably well. 

As for the excitation function, we have found that the resonance peak is reproduced 
within the errors without changing the isobar parameters in passing from pp pionpro duct ion 
to the pd one. We have also found that the magnitude of the curve {i.e. the normalization 
of the cross section) is sensitive to the nucleon-nucleon interaction used as input for the 3N 
bound state wavefunction, with a 25% difierence between the Bonn potentials and the Paris 
one. This sensitivity is larger by a factor of 5 with respect to the pp —> 7^~^d process.The 
results with the Paris potential were always smaller than those obtained with the Bonn 
potentials. 

We have then analyzed the angular dependence of the diEerential cross-section in the 
energy region spread around the resonance peak. At the peak there are practically no 










differences in the angular distribution with respect to the model selected for the nucleon- 
nucleon interaction. In moving away from the peak, the angular distributions begin to show 
some interaction dependence: the differences are larger at higher energies and backward 
angles. The calculation with the Paris potential seems slightly in better agreement with the 
general trend of the data at backward angles. Dehnitive conclusions, however, should be 
drawn only after the ttN interactions in s-wave will be also included. 

At the resonance peak we have singled out the role of non isobaric ttN interaction in 
p-wave (the lA term), analyzed the convergence of the calculation with respect to the 
number of intermediate three-baryon partial waves included, and considered the effect of 
the three-body dynamics in the nucleon-deuteron channel. It turns out that the size of 
these effects are all comparable. Indeed, the contribution of the lA term with respect to 
the total cross section ranges from 2 to 4 % depending on the model interaction which has 
been used. Differences of the same size are found when the calculation performed with a 
strictly necessary set of partial waves is compared with converged results. Finally, by means 
of the Faddeev-AGS formalism, it was possible to ascribe a 4 % effect to the contribution 
due to the three-nucleon dynamics in the nucleon-deuteron channel. It is obvious that 
this is equivalent to say that for the inverse pionabsorption reaction, hnal state interactions 
contribute with the same 4% amount. 

In conclusion, we have here conhrmed the expectations, previously formulated in Refs. 


r7| , |6^ , that a careful embedding of the basic pion absorption/production matrix elements 


in a Faddeev-based treatment of the few-nucleon dynamics is a very important tool for 
understanding the hadronic processes in nuclei at intermediate energies. Once this point 
has been settled, it is possible to move further on with the same approach and tackle other, 
more rehned experimental data underlining the pion few-nucleon systems. 

Such aspects are, e.g., the pionic absorption on the diproton |^, the wealth of exper¬ 
imental data involving polarization phenomena [|6|, |60| , |63| , |66| , the pion-induced reaction at 
energies around the pionic threshold |^, and, hnally, the meson-absorption coincidence 
experiments (vr’'',pp) at non conjugated angles with the connected puzzle of the “genuine” 
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three-body mechanisms [^. However, before that all these aspects can be theoretically 
disentangled, an important improvement is needed in our treatment (as well as in any other 
approach). This is connected with the role played by the pion-nucleon s-wave interaction 
(in both its isoscalar and isovector components) in multiple rescattering processes. Stated in 
other words, it is the role played by pion-nucleus final state interaction (or initial, depending 
on the selected direction in time). This aspect is still missing in the present treatment of 
the d(p,7r+)t process, but a flavour of its importance can be immediately percieved in the 
simpler case of the p(p,7r+)d reaction, by glancing at the dashed and dotted-dashed curves 
passing through the data in Figure 1. 
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FIGURES 


FIG. 1. Excitation function for d production (in microbarn) from pp collisions. The 
parameter p corresponds to the pion momentum (c.m.) divided by the pion mass. The full 
line is the result obtained with Bonn B potential and with p-wave pion-nucleon interactions 
which includes both isobar and non-isobar degrees of freedom. Practically indistinguishable 
from the full line, there is a dotted line which correspond to a similar calculation but with 
Paris potential. The other two curves (dashed, and dashed-dotted, respectively) correspond 
to the inclusion of pion-nucleon initial state interaction in s-wave and differ between each 
other for the NN potential being used. The Bonn B in the dashed case, Paris in the dotted- 
dashed case. 

FIG. 2. Excitation function of pion production via proton-deuterium collisions (in mi¬ 
crobarn). The calculation includes also A excitation and has been performed with the same 
parameters as in Figure 1. 

FIG. 3. Differental cross-section for tt’*' production in proton-deuterium collision. 

FIG. 4. TT■'■-production differential cross-section at p = 1.36. The full line shows the re¬ 
sults with inclusion of nucleon-nucleon initial-state interactions. The dotted line represents 
the plane-wave calculation. The dashed line shows the results with a limited number of 
channels. The dashed-dotted line shows the contribution of the sole non-isobar degrees of 
freedom, magnihed by a factor 10. All calculations were performed with the Paris potential. 
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FIG. 5. Unpolarized differential cross-section for 7r° production in proton-deuterium 
collisions. Circles and triangles represent measurments taken at Saturne at forward 
and backward (c.m.) angles, respectively. The solid and broken curves are forward-angle 
calculations with Bonn B and Paris potentials, respectively, the dashed and dotted-dashed 
curves are the corresponding backward-angle results. 
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TABLES 
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TABLE I. Example of three-nucleon partial waves included in our calculation. The notation 
is A/, where s, I, are j are spin, orbital, and total momentum of the pair, while A and I 

are orbital and total momentum for the spectator nucleon. This set of 82 states corresponds to 
the inclusion of 4 two-nucleon states. In our actual calculation 464 three-body partial waves have 
been considered, corresponding to the inclusion of 18 two-nucleon states. 
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Bonn A 

Bonn B 

Paris 

lA 

0.58 

0.72 

0.93 

PWa 

30.1 

28.0 

21.7 

PWb 

30.7 

28.7 

22.5 

ISIa 

31.3 

29.0 

22.5 

ISIb 

32.0 

29.9 

23.4 


TABLE II. Calculated production cross sections (in microbarns) at ?? = 1.36 for various NN 
potentials. The Impulse Approximation (first row) represents a bare plane-wave calculation with¬ 
out A rescatterings. The second and third rows include also the A degrees of freedom but use a 
plane-wave approximation in the incoming channel. They differ between each other for the number 
of partial waves included in the intermediate states. In case a, 4 two-nucleon states have been 
included and only S'-wave orbital momentum for the intermediate A-nucleon state has been con¬ 
sidered. In case b, 18 two-nucleon states and intermediate A-nucleon states in S, P, and D waves 
have been taken into account. In the last two rows, where the same number of partial waves have 
been included as in the previous two, the p-d interactions in the initial state have been included 
via a Faddeev calculation. 
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partial-wave 

PEST 

BAEST 

BBEST 
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4 
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CO 
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4 

3P2 - ^F2 

7 

5 

5 

TABLE III. Ranks of the two-body partial-waves of the Paris, 

in the EST-representation. 

Bonn A, and Bonn B potentials 

BAEST 


BBEST 

PEST 

-8.284 


-8.088 

-7.3688 


TABLE IV. Calculated binding energies for the triton. The total angular momentum of the 
two-body subsystem was restricted to j < 2. This led to 18 three-body channels in the chan¬ 
nel-spin coupling. 
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Fig. 1 Canton PRC 
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